Deviations from the Gaussian distribution of mesoscopic conductance fluctuations 
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The conductance distribution of metallic mesoscopic systems is considered. The variance of this 
distribution describes the universal conductance fluctuations, yielding a Gaussian distribution of the 
conductance. We calculate diagrammatically the third cumulant of this distribution, the leading 
deviation from the Gaussian. We confirm random matrix theory calculations that the leading 
contribution in quasi-one dimension vanishes. However, in quasi two dimensions the third cumulant 
is negative, whereas in three dimensions it is positive. 
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I. INTRODUCTION 

About, rten years ago the universal conductance fluctuations (UCF) were discovered in mesoscopic, metallic 
samplesErEl. The electronic conductance of these samples shows reproducible sample to sample fluctuations. The 
fluctuations are called universal because their magnitude is independent of the sample parameters such as the mean 
free path £, and the average conductance (g). The dependence on sample dimension is weak. Studies have mainly 
focused on the variance of the fluctuations, which is the leading and the universal part of the fluctuations. 

The conductance being a random variable showing such large fluctuations, one realized that one should consider 
its full distribution. It was soon clear that the first higher cumulants of the conductance are proportional tea 

(g») c ex (g) 2 -", n < go, (g) » 1. (1) 

Here the conductance is measured in units of e 2 /h, go is the mean conductance at the scale £, (. . .) denotes the ensemble 
averaging, and the subscript c indicates cumulants. In the metallic regime far from localization, where (g) 3> 1, the 
higher cumulants are thus small, and the distribution of the conductance is therefore roughly Gaussian. However, 
for n > go the decrease in magnitude of cumulants as described b^ Eq. ([!]) is changed into a very rapid increase 
(oc expfg^n 2 ]). This leads to the log-normal tails of the distributionQ. With increasing disorder, the log-normal tails 
become more important. Although the calculation of the full conductance distribution on the threshold of localization 
((g) ~ 1) is today out of reach, it is quite plausible that the whole distribution crosses over to a log-normal shape in 
the strongly localized regime (see Refs. || and ^ for a discussion). Indeed, it is well known that in the strongly localized 
regime in one-dimension the conductance is given by the product of transmission amplitudes, yielding a log- normal 
distributiontm. 

Although high order cumulants govern the tails of the distribution, they do not affect it near the center. A deviation 
from the Gaussian distribution near the center should be revealed, first of all, in the lowest nontrivial cumulants. An 
important step in this direction was the recent calculation of the third cumulant of the distribution using random 
matrix theory by MacedcO. He found for the orthogonal (fi — 1) and symplectic ensemble (fi = 4) that the third 
cumulant of the conductance is proportional to 1/g 2 , thus the leading term given by Eq. (1) vanishes (see for instance 
Ref. |lj for the definitions of the ensembles). For the unitary ensemble (f3 = 2) even this sub-leading term vanishes. 
The physical reason behind this is not clear. (We recently learned that the same result in quasi one dimension was 
found by TartakovskiO by the scaling method described in Ref. |lj.) However, random matrix theory is only valid in 
quasi one dimensional systems. Therefore, it was not known whether this cancelation is unique to one dimension, or 
holds also in higher dimensions, which might indicate an overlooked symmetry of the system. That such a possibility 
exists was also suggested by the fact that the leading-order contribution to the third cumulant of the density of states 
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vanishes in two dimensions!! The question whether or not the cancelation holds in two and three dimensions was a 
major motivation for this work. ■— ■ 

Let us finally mention that related to the conductance distribution, recently third cumulantst3 and the full dis- 
tribution functionsEj of two related transmission quantities were calculated. These quantities, the speckle transmis- 
sion and the total transmission, were measured in scattering experiments with electro-magnetic waves (light and 
microwaves)E3ii3. In the regime considered, these distributions are independent of dimension. On a diagrammatic 
level thiSpCan be understood: The diagrams for the cumulants of the speckle and the total transmission are loop-less 
diagramsE£l. The diagrams for the cumulants of the conductance, however, contain loops, and will therefore be, a 
priori, dimension dependent. For the electro-magnetic waves the Laadauer approach was used, however, for conduc- 
tance properties this yields the same results as the Kubo approachEjO Therefore, the calculations presented below 
are also valid for classical waves. 

In the present paper we will calculate the third cumulant of the conductance for any rectangular geometry in the 
metallic regime in quasi one, quasi t wo, and three dimensions. In section [n] we review the formalism for ca lculat ing 
conductance fluctuations. In section III it is described how the leading diagrams can be found. In section [II A we 
calculate the third cumulant and discuss its dimension dependence. After considering the effects of inelastic scattering, 
we end with a conclusion. 



II. CONDUCTANCE FLUCTUATIONS 



In this sectiop-we introduce our diagrammatic approach, to a large extent following the detailed paper of Lee, Stone, 
and FukuyamaEl We consider a sample with a rectangular geometry with sizes L x x L y x L z ; the conductance is 
measured across the z-direction. The scatterers are isotropic point scatterers. The scattering is calculated in second 
order Born approximation, which assumes weak scattering. (The present calculation can be easily-,,e?j:tended beyond 
second order Born, this will yield the same results provided the correct mean free path is takenr^r 3 !.) We consider 
the metallic, mesoscopic regime, which is characterized by the inequalities 1/k <C I <C L z , where k denotes the Fermi 
wavenumber; £ is the mean free path, related to the scatterer density n and scatterer strength v as I = 47r/n\v\ 2 . The 
relaxation time is given by r = I /(2k). Incoherent scattering is first neglected, that is, the incoherence length is much 
larger than any system size, L m 3> L Xl L y , L z . 

The conductance and its cumulants are calculated with the Kubo formula. The transport through the sample is 
given by the so called ladder diagrams or diffusons, describing the multiple scattering of the electrons on the scatterers. 
In the conductance diagram the field creates an electron-hole pair at a current vertex somewhere in the sample, after 
some diffuse propagation this pair annihilates. This diagram thus contains a bubble, the corresponding propagator C 
obeys the diffusion equation: 

-V 1 £(r 1 ,r 2 ) = ^-l^«5(r 1 -r 2 ), (2) 

where D is the diffusion constant. (Correspondence with units often used with classical waves is found by identifying: 
K = l,m = 1/2, v = fc/(47r 2 ), D = 2H/3.) 

As mentioned we restrict ourselves to the conductance in the z-direction. The current is restricted to the z-direction 
by imposing fully reflecting boundaries in the x and y direction, and fully conducting boundaries in the z-direction: 

ru \\ n d£(ri,r 2 ) dC(r 1 ,r 2 ) 

M r l; r 2jUi=0,zi=L 2 — U, — \ Xl =0,xi=L a — g— |yi=Q,« 1 =i 1) - U, W 

and similarly for the r 2 dependence. For rectangular geometries it is useful to write the solution of the diffuson 
propagator in (discrete) momentum space. Following Ref. |2l| (in present work there is an extra factor 1/(2ttut) in 
the diffuson), we write the diffuson as 

oo oo oo 



f( r r )- 1 V V V" Q'( r l)Q'( r 2) (a) 

£(ri ' ra) - 2Wr7 2. 2.2. ^ (4) 

i s — 1 ix— i y — 

where i is the momentum vector (ix,iy,iz)- The Qi's are the normalized, orthogonal cigenfunctions 
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Note that these eigenf unctions are not properly normalized if i x = or i y — 0, in which case one has to replace 
cos (t^) or cos ( Z x7) by l/y/2. The eigenvalues A t in Eq.(@) are 



(0) 



The longitudinal momentum takes positive integer values; the transverse momenta i x and i y can in addition also 
be zero. In quasi two dimensional (2D) geometries one has L x -C L z and only the terms with i x = contribute in 
the eigenvalues. In quasi one dimensional (ID) samples also L y -C L z , and i y is essentially restricted to zero. On the 
other hand, for very wide geometries i x (L z /L x ) 2 becomes a continuous variable and the sum over the x- variable (or 
y- variable) becomes an integral. The average value of the dimension-less conductance (g) is 

(g) = 2wD^. (7) 

The inverse of this, l/{g), will be a small parameter in our diagrammatic expansion. 

A. Hikami-boxes 

For the calculation of the second and higher order cumulants p£ @~ we need the vertices describing the interference 
between two diffusons. These vertices are known as Hikami-boxesE-lEj. Formally they arise from the spatial derivatives 
d/dr a and d/dr'p in the Kubo formula for cr Q ^(r,r'), where the indices a, (3 — x,y,z label the directions of the 
incoming and outgoing current. We have drawn the boxes in Fig. |l|, the corresponding expressions are labeled H a to 
Hf, respectively: 

H a = 4n^ S -^f, (J 

H b = 2W^, (8b 
k 

H cl = 4niisT 4 q%-, (8c 
3 

H c2 = -H cl , (8d 

H d = 2ttvt a D [-2qi • q 3 - 2q 2 • q 4 - (qi + qa) • (<J2 + q4)] , (8e 

He = Attvt 3 ^, (8f 
3 

H f = -2itvt 5 ^. (i 

The q's are the momenta of the diffusons (the wavy lines), as numbered according Fig. [l| As the conductance is 
calculated in the z-direction, a and j3 are taken z. Likewise, only the z-component of the momentum in H c will be 
taken into account. H c comes in two flavors: H c \ and H C 2- H C 2 has an additional minus sign due to the reversed 
orientation of the advanced and retarded propagators of the current vertex. We used that the diffusons vary slowly 
on length scales of the mean free path, i.e. Finally, other boxes are sub-leading. 



B. UCF calculation 



For clarity we briefly present the UCF calculation using this formalism. The UCF-diagrams are connected diagrams 
containing two conductivity bubbles and thus four current vertices. Its diagrams are shown in Fig. |. The UCF 
diagrams contain two four-boxes (H a or H b ) and two diffusons. For the orthogonal ensemble the combinatorial 
factor for the upper diagram is 4. The lower diagram has combinatorial factor 2, and the two Hikami vertices bring, 
comparing Eq. (|8a|) to Eq. (Sb), an additional factor of 2 2 . Thus one finds a pre- factor 12 for the sum of all diagrams 
with respect to a single diagram of the upper topology. (In the unitary ensemble cooperons do not contribute and 
this would be a factor 6.) For each incoming or outgoing current vertex there is a factor 2k/ L z . It arises from the 
expression l/(mL z ) — > 2/L z , that can be read off from Eq.(2.5) in Ref. [2l], and the wavenumber k that originates from 
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each derivative in the Kubo formula. Next, we use the Fourier decomposition of the diffusons Eq.(Q) and interchange 
the sum over the eigenvalues with the spatial integrals. The two diffusons interfering at the Hikami-box yield the 
orthogonality relations: J drQi(r)Qj(r) = s\j. One finds 



h3 



7T* — V" L l L l 



In quasi one dimension this yields ^| 5Di°=i *z 4 = T51 which is the well-known result for the UCF in quasi one 
dimension0. 



III. THE THIRD CUMULANT 

We now study the third cumulant. First, we rewrite Eq.(0) into a relation for the relative cumulants: <jy£ ^ 

(g) 2 ~ 2n . The inverse power of (g) on the right hand side can be interpreted as the number of Hikami four-boxes 
in the diagrams. Thus the third cumulant diagram contains four Hikami four-boxes. Indeed, it proves impossible 
to create {g 3 ) c diagrams with less boxes. Diagrams with more boxes are sub-leading as they are of higher power of 
1/(17), which is a small parameter. In diagrammatic approaches it is often tricky to find all the relevant diagrams; 
already for the much simpler set of UCF diagrams there was considerable discussion in the literature. We used the 
following considerations to find the set of diagrams for the third cumulant: The diagrams have six current vertices: 
three incoming ones and three outgoing ones. The diagrams contain four 4-point vertices and each vertex has at least 
two diffusons attached to it, see Fig. |[ This leaves two possibilities: there are three boxes with each two diffusons and 
one box with four diffusons; or, are two boxes with two diffusons and two boxes with three diffusons each. It is now an 
easy exercise to see that there are only four possible basic topologies. They are drawn in Fig. |[ The dots represent 
the 4-point vertices; the lines represent diffusons. The next step is to insert the Hikami-boxes and connect them in 
all possible ways to the diffusons. This yields many diagrams, yet not in all diagrams the outgoing electron-hole pair 
is the same as the incoming electron-hole pairing. Such diagrams are not the product of three conductance bubbles; 
they do not contribute to the (g 3 ) c process and are left out. We end up with the diagrams represented in Fig. |^. 
In the figure there are also diagrams of the same order with one six-box and two four-boxes. They can be obtained 
by contracting one diffuson in the diagrams with four-boxes. In physical terms the six-box diagrams correspond to 
processes where after an interference process, the amplitudes do not combine into a diffuson. Instead, the amplitudes 
interact again directly without being scattered. 



A. Evaluation of the diagrams 

We calculate the diagrams for quasi-lD, quasi-2D, and 3D geometries. The calculation goes similar to the UCF- 
calculation, that is, the sums over the momenta are taken outside the spatial integrals. If two internal diffusons 
interfere, the spatial integrals bring the orthogonality relations as used for the UCF calculation. However, two 
complications arise. First, at the Hikami vertices not only two, also three, or four diffusons interfere (as is seen from 
from the Figs. |^ and [IJ), this gives more complicated structures. Secondly, Hikami-boxes H c and Hd are g-dependent, 
corresponding to spatial derivatives of the diffusons. We therefore introduce the matrices R, S, and T, which describe 
the interaction of the diffusons at the vertices. 

In the diagrams e.i), e.ii), e.v), and e.vi) of Fig. ^ four diffusons, two with momentum i, and two with momentum 
j interfere. The x-components of the diffusons couple as 



1 + 2(1-^,0)^ 



(10) 
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The y-components yield exactly the same integral, see Eq.(j|), and also the z-components yield the same expression. 
These diagrams will thus be proportional to Ri z j z R^ j x Ri y j y . 

In most other diagrams three diffusons interfere at a Hikami box H c \, this box brings a derivative d/d z . The matrix 
S describes the effect of this derivative on the z-component, the T describes the integrals over x and y components. 



3/2 r L z 



L 



dz 



sin 



j z zn 



sm 



L z 



/2L 



,j* + k z + Sj z A z +k z + Sk z A z +j z ) 



2x ,ix 



2 

L x 
1 



3/2 



dx cos 



IzZ-K 

~T7 



cos 



cos 



L x 
1 



(^,0^,fe s +^,o^,fe„+^,o^^ -2^,0^,0^,0) (if ix, jx, or k x = 0) 



: (^*x,jx+fcx ")"^jx,ix+fex +^fex,4x+jx) 



(else) 



(11) 

(12) 
(13) 



The box H C 2 has a minus sign resulting in —Si zj j z _k z T t 



^x tjx i^x iy ijy 



The diagrams can now be written in terms of the R, S, T. The pre-factors, the combinatorial factors, and the sum 
over the momenta are included. The third cumulant is given by the sum of diagrams (g 3 ) c = where the F's 

are, 



P . 4. P . . = _ 1 AA r \ " 1 <? 2 . , T 2 , T 2 

a - 1 r a.%i \3\ \ a » i ,^,fc i - t «x,Jx,fex- t » v ,. 



(14a) 



&.» 



i 7 ? 



Ll V2 



L *a: iJx iy tjy 



It 



2 

E T2T2 



^f2 ^2 /^>2 

^2^2^ i>z jj z i^z "ixjjxj^x iyijyjky 

i,j,k 1 i 



Rc.ii ~t~ -Fcili 



120r^ \2\2\ Sk x ,u ,jJ^i x ,j x ,kJ^iy 



I A A " A * 

i.j,k L j 



jy i^y 



F, 



d.i 



F R 



Fd 



Fd 



Fa An] 



96r E \2\2\ Sizj*>kzSj z ,i z< kzT? xt j xtkx T? yt: j ytky 
864r E ylYIT s <.,i..*. 5 i.^,*. T <-,<..fc- T i-J-,*- T *,^,^ T i»A.*, 



. A 2 A 2 A, 



Fe.ii H~ F ev ~\~ Fe.vi ~ 367*^^ \2\2 ^^zjz^ixd 



A A * 



(14b) 
(14c) 

(14d) 

(14e) 

(14f) 

(14g) 



The sums runs over all allowed momenta indices. The factor r contains the pre-factors of the diffusons and boxes and 
a factor (2k/L z ) 6 for the incoming and outgoing current vertices, 



2ttv DL* 



(15) 



Finally, the sums still contain a factor L 7 z /(L x L y ), yielding multiplied by r a pre-factor (g)^ 1 . So indeed the third 
cumulant is proportional to the inverse dimension-less conductance as predicted. The choice of ensemble is reflected 
in the combinatorial factor. For simplicity the combinatorial factors were calculated for the unitary ensemble, finally, 
for the orthogonal ensemble all pre-factors will be four times larger. The e.iii) and e.iv) diagrams are divergent in two 
and three dimensions, as the sum Yl^ ^ ^ XJ 1 is logarithmically divergent in 2D and linearly divergent in 3D. This 
divergence is, however, exactly canceled by a similar divergence in the b.i) and b.ii) diagrams. Eq. (14b) presents the 
combined, finite expression. 
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It is interesting that due to the finite system size momentum conservation apparently does not hold, for example, 
the momentum of the middle diffuson in the d) diagrams of Fig. |J need not be zero. This is a result of the mirror 
terms in the diffusons, which are present to fulfill to the boundary conditions. Though such terms would be absent 
in the bulk, they never vanish in our geometry. 

In the quasi one-dimensional case the diffusons are simple linear functions of the z— coordinate 



1 



\(z 1 ,z 2 )[L z - max(zi,z 2 )] 



2kvt 2 D 



(16) 



allowing an analytically treatment. The diagrams can now even be calculated directly without introducing the 
momentum representation at all. As a check we also calculated the diagrams using this representation. We find, using 
cither method, for the sum of the diagrams 







(quasi ID). 



(17) 



Thus the leading contribution to the third cumulant in one dimension vanishes. This confirms the random matrix 
theory resultliil diagrammatically. 

In higher dimensions the sums were performed numerically. In the numerical evaluation the sums over three sets of 
momenta are implicitly reduced to sums over two sets. This is a direct consequence of the fact that there are only two 
independent momenta present. The numerical evaluation remains, however, quite involved as in three dimensions one 
still has to sum over six variables and the convergence is quite slow. We find in the orthogonal ensemble for square, 
and cubic samples, resp. 



<.9 3 > 
(9 3 ) 



-0.0020( 5 ) _1 
-0.0076(s) _1 



(quasi 2D) 
(3D) 



(18) 



The results for rectangular samples are given in the figures ^| and ^|, where we multiplied the third cumulant by the 
average of the dimension-less conductance. The third cumulant for wide slabs (L x 3> L z , L y L z ) is proportional 
to (L x L y /(g)). In the figures due to the multiplication by (g) there is proportionality to (L x L y ) 2 / L\ for wide slabs. 
For very narrow slabs one sees that the correct quasi 2D and quasi ID limits are recovered. Note that third cumulant 
passes through zero when going from 2D to 3D, this happens if the sample has the size 0A6L Z x L z xi 2 . 

B. Inelastic scattering 

Incoherent scattering was neglected in above calculations. In realistic system, however, incoherent scattering can 
be present. This is especially true at non-zero temperatures where electton-phonon interactions will occur. This 
mechanism was included in the description of Lee, Stone, and FukuyamaEil. The inverse inelastic scattering time 
induces a positive shift of the diffuson eigenvalues 



K ~ LI 



TT 2 L, 



(19) 



The case that the incoherence length is much smaller than the system sizes is of particular interest. The effect on 
the third cumulant can now be estimated simply by considering the sample as being made up of small samples of 
dimension Lf n with independent conductance distributions. The conductance of such a sample is in 3D 



^rc/^in L y j Li] 

E E 



-go 



L,/L ia 

E (sor 1 



L X Ly 



E^o, 



(20) 



where g denotes the conductance of a individual coherent cube. The relative cumulants thus scale as 

(9 n )c_f Lf n \ n (g8) c 



{g) n \L x L y Lj ( 5o )"' 
The relative third cumulant reduces by a factor Lf u /(L x L y L z ) 2 . As a result the combination 



(21) 



G 



(g 3 )c(g) , 99 x 

is independent of L m , so that this quantity is universal. As the incoherence length is the same in all directions, the 
coherent parts are essentially cubic. Therefore, the quantity will roughly tend to its value for coherent, cubic samples. 
The precise value can be obtained by properly including the incoherence effects in the calculation as indicated above. 
Because of this universality, the quantity is useful experimentally. We even expect the full conductance distribution 
to be universal in this regime, that is, independent of geometry and incoherence length. 



IV. CONCLUSION 



We have considered a mesoscopic sample in the metallic regime and we calculated the third cumulant of the 
conductance distribution. Naive scaling predicts that the third cumulant should be proportional to I /(g). In quasi 
one dimension we confirm the absence of this leading contribution, as was found by MacedcO. In two and three 
dimensions, however, this is cancelation is not present; the leading contribution to the third cumulant is negative 
in two dimensions and positive in three dimensions. The fact that the third cumulant changes sign is surprising. 
The third cumulapt is also known as the skewness of a distribution. In analogy with the third cumulant of the 
total transmission!^ or if the distribution would be tending to log-normal, one would have expected a positive third 
cumulant of the conductance. Instead, we find that all possible values occur: negative, positive and zero. We have 
no explanation for this. 

To the best of our knowledge there exists no experimental work where the conductance distribution is discussed. 
Experiments, cither electronical, or using classical waves, or numerical simulations, could enlighten present results. 
Predicted values are small, but should detectable in electronic systems with moderate values of (g). Also observation 
of the conductivity distribution as a whole would be very interesting. 
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Figure 1 

FIG. 1. Hikami four-point and six-point vertices used. The wavy lines denote diffusons or cooperons. The current vertices 
are at the unoccupied corners of the polygons and are labeled by a and p. 




FIG. 2. The UCF or {g 2 ) c diagrams. 




FIG. 3. All basic topologies for the third cumulant diagrams. The wavy lines are diffusons or cooperons. The dots represent 
Hikami four point vertices. Detailed inspection of these diagrams will yield the diagrams of Fig. W. 
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FIG. 4. The calculated set of leading diagrams for the third cumulant of the conductance. The diagrams are derived from 
Fig. ^; one sees that structure is the similar. The diagrams e) are obtained by contraction of diffusons in the diagrams a), b), 
c) and d). 

Third cumulant of conductance 

Width dependence quasi 2D 
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width ratio: LJL z 

FIG. 5. The third cumulant of the conductance multiplied by the average as a function of the transversal size in a 2D sample. 
Although the leading term vanishes in quasi ID, it is non-trivial in 2D geometries. 
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FIG. 6. The third cumulant multiplied by the average conductance in 3D sample plotted against the transversal size L x . 
The geometry is L x x L z x L z for the solid line. Note that, the sign changes in going from a quasi-2D to a 3D sample. The 
dashed line gives third cumulant in sample with geometry L x X L x x L z . 
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